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w  ON  THE  EFFECT  OF  SUBSONIC  TRAILING  EDGES  ON  DAMPING  IN 

o 

ROLL  AND  PITCH  OF  THIN  SWEPTBACK  WINGS 
IN  A  SUPERSONIC  STREAM 
By  Herbert  S.  Ribner 


SUMMARY 

The  principal  effect  of  subsonic  trailing  edges  on  the  damping 
in  roll  and  pitch  of  thin  sweptback  wings  in  a  supersonic  stream  is 
evaluated  with  the  aid  of  some  conical  and  quasi -conical  flows  pre¬ 
viously  derived.  This  effect  is  expressed  in  the  form  of  approximate 
correction  terms  to  be  added  to  the  corresponding  expressions  that 
are  obtained  when  the  trailing-edge  disturbance  is  ignored.  The 
results  are  limited  to  those  plan  forms  and  Mach  numbers  for  which 
the  trailing-edge  disturbance  does  not  extend  beyond  the  leading  edge 
and  for  which  the  area  of  mutual  interference  between  tips  and  trail¬ 
ing  edge  is  not  large.  Practical  applicability  is  subject  to  the 
limitations  of  linearized  potential  theory. 


INTRODUCTION 

At  supersonic  speeds  the  flow  over  a  sweptback  wing  is  closely 
related  to  the  flow  over  a  triangular  or  delta  wing  (fig.  1(a))  with 
the  same  leading-edge  sweep  angle.  Thus,  in  figure  1(b),  the  flow 
characteristic  of  the  delta  persists  out  to  the  tip  Mach  cones. 

Within  the  tip  Mach  cones  the  flow  is  altered.  (See  reference  1.) 

At  lower  supersonic  speeds  (fig.  1(c)),  the  conical  Mach  wave  from 
the  apex  of  the  trailing  edge  encloses  a  portion  of  the  trailing  edge. 
Here  also  the  basic  delta  flow  is  altered.  (At  these  speeds  the  com¬ 
ponent  stream  velocity  normal  to  the  trailing  edges  is  subsonic  and 
the  edges  are  termed  "subsonic.") 

The  altered  flow  in  the  tip  and  trailing-edge  regions  may  be 
considered  to  result  from  superposing  on  the  basic  delta -wing  flow 
certain  necessary  additional  flows.  These  additional  flows  have  the 
function  of  canceling  the  delta -wing  lift  outside  the  wing  tips  and 
behind  the  trailing  edge.  These  cancellation  flows  yield  the  dis¬ 
turbance  pressures  in  the  tip  regions  and,  if  the  edge  is  subsonic, 
in  the  trailing-edge  region. 
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According  to  figure  1(c) ,  the  tip  and  trailing-edge  disturb- 
ance  regions  overlap.  In  the  cross-hatched  upper  portion  of  the 
overlap,  for  each  tip,  the  trailing-edge  disturbance  and  the  tip 
disturbance  merely  superpose.  Secondary  disturbances  in  the  black 
lower  portion  of  the  overlap  sure  caused  by  mutual  interference 
between  the  tip  and  the  trailing  edge.  Under  fairly  general  cir¬ 
cumstances,  these  secondary  disturbances  have  been  found  (reference  1) 
to  contribute  a  relatively  unimportant  amount  to  the  integrated  lift 
and  moment.  Accordingly,  the  secondary  disturbances  are  neglected 
herein  and  the  tip  and  trailing-edge  disturbances  are  assumed  to  be 
independent  and  to  superpose  within  the  entire  region  of  overlap. 

The  problem  of  the  tip  disturbance  has  received  a  rigorous  treat¬ 
ment  in  references  1  and  2  and  a  simplified  approximate  treatment  in 
reference  3.  The  discussion  herein  is  limited  to  the  trailing-edge 
disturbance. 

The  complete  flow  necessary  to  cancel  the  lift  behind  a  sub¬ 
sonic  trailing  edge  is  difficult  to  evaluate.  It  is  shown  in  refer¬ 
ence  1,  however,  for  the  case  of  angle  of  attack,  that  the  major 
part  of  the  disturbance  pressure  in  the  subsonic -trailing-edge  region 
is  obtained  with  only  partial  cancellation.  (The  complete  cancel¬ 
lation  is,  however,  carried  out  to  close  approximation  in  refer¬ 
ence  1.)  This  partial  cancellation  is  afforded  by  a  certain  conical 
flow  (designated  flow  I)  in  the  manner  shown  in  figure  2(a).  The 
flat  portion  of  the  load  distribution  in  section  A-A  represents  a 
uniform  negative  lift  in  the  triangular  region  behind  the  trailing 
edge;  its  magnitude  is  chosen  to  equal  the  value  of  the  basic  delta - 
wing  flow  at  the  midspan.  The  outer  portions  constitute  the  sought - 
for  pressure  disturbance  on  the  wing  in  the  trailing-edge  region. 

Similar  concepts  can  be  employed  to  evaluate  in  simple  fashion 
the  principal  portion  of  the  subsonic -trailing-edge  pressure  dis¬ 
turbance  in  the  case  of  rolling  or  pitching.  For  rolling  motion  the 
basic  lift  to  be  canceled  is  antisymmetric.  The  appropriate  partial 
cancellation  is  afforded  by  a  certain  flow  III  in  the  manner  shown 
in  figure  2(b).  For  pitching  motion  the  basic  lift  to  be  canceled 
increases  linearly  downstream.  The  appropriate  partial  cancellation 
is  afforded  by  superposition  of  a  constant  negative  lift  (flow  I) 
and  another  lift  (flow  IV)  that  increases  in  proportion  to  the  dis¬ 
tance  downstream  of  the  apex  of  the  trailing  edge  (fig.  2(c)). 

The  detailed  derivation  of  these  several  cancellation  flows  or 
"wake  corrections"  and  others  is  carried  through  in  reference  4.  The 
present  report  is  concerned  with  the  application  of  some  of  these 
results  to  the  determination  of  damping  in  roll  and  pitch  for  thin 
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swept  wings  with  subsonic  trailing  edges.  These  applications  +»Tp» 

‘  the  form  of  correction  terms  to  be  added  to  the  corresponding  expres¬ 

sions  that  are  obtained  when  the  trailing-edge  disturbance  is  ignored. 

H 

o  The  effect  of  subsonic,  trailing  edges  on  the  ^wmp-tng  in  roll  of 

sweptback  wings,  which  is  part  of  the  present  subject,  has  also  been 
treated  recently  in  reference  2.  The  treatment  therein  achieves 
nearly  exact  cancellation  of  the  lift  behind  the  trailing  edge  by  the 
superposition  of  infinitely  many  infinitesimal  flows.  The  procedure 
is  a  development  of  the  approach  used  in  reference  1.  The  present 
treatment  differs  in  employing  approximate  cancellation  by  super¬ 
position  of  a  single  flow.  This  simplification  permits  giving  the 
over -all  result  compactly  in  closed  form. 

The  results  presented  herein  are  applicable  for  those  super¬ 
sonic  speeds  for  which  the  trailing-edge  disturbance  does  not  extend 
beyond  the  leading  edge.  At  the  lower  Mach  numbers  for  which  the 
leading  edge  is  enveloped  further  correction  is  required.  This  low 
Mach  number  problem  is  studied  for  the  case  of  angle  of  attack  in 
reference  5. 


CORRECTION  FOR  DAMPING  IN  ROLL 

For  a  sweptback  wing  with  subsonic  leading  edges  rolling  with 
angular  velocity  p,  the  basic  delta-wing  lift  distribution  to 
be  canceled  behind  the  trailing  edge  is  (equation  (26),  reference  6) 

2 

2(x+cr)ypm  l(m) 

°P,a - —r  :  (l) 

W  \j  (l+or)2«2-S2y2 


where 

I(m)  =  - Hl-m) - 

( 2-m^  )E ' (m)  -m^K'  (m) 

and  x  is  measured  from  the  apex  of  the  trailing  edge.  (All  symbols 
are  defined  in  appendix  A.)  This  load  distribution  has  the  anti¬ 
symmetric  character  shown  in  figure  2(b).  For  the  inexact  cancel¬ 
lation  employed  herein,  the  distribution  is  approximated  by  its 
tangent  at  y=0: 
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y 


S  2Pml(m2  y  (2) 

pv 

(The  limitations  on  the  use  of  this  approximation  are  considered  in 
the  DISCUSSION.) 

The  load -distribution  equation  (2)  can  be  canceled  by  flow  III 
of  reference  4,  which  can  be  written  as 

Cp  =  Cy  (3) 


t 


where 

C  =  constant 

in  the  region  0  <.  |  o|  S.  n  behind  the  trailing  edge  of  the  wing. 
Thus  the  choice 

2pml(m)  /.v 


achieves  the  desired  cancellation. 

Cancellation  flow  III  includes  also  a  lift  disturbance  ahead  of 
the  trailing  edge  in  the  region  n  <.  |al  <1.  Deference  4  gives  for 
this  disturbance,  in  the  present  notation, 


o 

1  -  \ 

f(i-a2)(a2-n2) 

3 

[jj(k)-] 

n2K(kj] 

J 

where 


C  =  same  constant  as  in  equation  (3) 


$  =  sin"-1 


(5) 
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k  =  *\j  1-n2 

C  =  Py/x 

The  subsonic -tralling-edge  correction  to  the  damping  in  roll  results 
from  integrations  to  obtain  the  contribution  of  equation  (5)  to  the 
rolling  moment  of  the  wing. 

The  correction  is  evaluated  as  follows:  Equation  (5)  is  written 
in  the  form 


Cp  «  Cxf(o) 

Then  the  incremental  rolling  moment  on  an  elementary  triangular  region 
(hatched  area,  fig.  3)  is  given  by 


/ 


Pb 
p2  a 

yCxf(a) 


x=0 


cp2b4f(o)  da 

6405 


« 


The  incremental  rolling  moment  experienced  by  the  finite 
triangular  region  n  <  O  <  1  of  figure  3  is  given  by 


AL* 


2  4 
CP  b 

64 


^n 


(6) 


The  function  f (a)  is  the  expression  in  braces  in  equation  (5). 
Thus  the  integrand  in  equation  (6)  involves  elliptic  integrals  and 
related  functions.  The  integration  may  be  effected  without  diffi¬ 
culty  with  the  aid  of  the  elliptic -function  substitutions  of 
appendix  B.  The  result  is 


*  l-P2 
4  E(k)-n2K(k)_ 
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The  following  standard  definitions  are  convenient: 

E(k)  =  E*(n) 

K(k)  =  K*  (n) 

With  these  substitutions,  equation  (7)  becomes 

AL*  =  _  Cpb4 
|  pY2  64n 


1  -  S 


1-n 


4  E’(n)-n2K'  (n) 


(8) 


The  value  of  C  given  by  equation  (4)  is  to  be  substituted  back 
in  equation  (8).  Then  the  expression  may  be  converted  to  coefficient 
form  by  division  by  Sb  and  by  pb/2V.  Multiplication  by  a  factor 
of  2  accounts  for  both  wing  panels  on  the  complete  wing.  The  final 
result  is  (with  both  sides  multiplied  by  P) 


A (PC,  ) 
ZP 


APml(m) 

8n 


1  -  - 


l-n‘ 


4  E’  (n)-n2K‘  (n) 


(9) 


in  which  the  leading-  and  trailing-edge  sweep  parameters  m  and  n, 
respectively,  are  related  by 


4 


11  4  1-X 

n  ”  m  ”  AP  1+X 


(10) 


where  A  is  the  aspect  ratio  and  X  the  taper  ratio  of  the  wing. 

Equation  (9)  gives  a  correction  for  the  damping -in -roll  deriv¬ 
ative  C7  when  the  trailing  edge  of  the  wing  is  subsonic.  This 

correction  is  algebraically  additive  to  the  value  of  PC,  computed, 

P 

as  in  reference  3  (in  effect),  by  integration  of  the  basic  delta¬ 
wing  lift  distribution  together  with  the  primary  tip  correction. 


Applicability  of  equation  (9)  is  limited  to  those  plan  forms  and 
supersonic  speeds  for  which  the  trailing-edge  disturbance  does  not 
envelope  the  leading  edge  of  the  wing.  This  limitation  is  a  conse¬ 
quence  of  the  shape  of  the  region  of  integration  assumed  in  fig¬ 
ure  3.  The  limitation  may  be  expressed  by 


1  <  -  <  1  + 
n 


±  X 
A|3  1+  X 


* 


(lla) 
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or  alternatively  by 


m 


AP(l+\) 


(lib) 


CORRECTION  FOR  DAMPING  IN  PITCH 

Pltc^oS8!tf*8~k  7*f  ^  8“l8°“la  leading  edge 

—  -  £2VS2 


PyA 


4qG(m)  [2(x+cr)2m 2~p2y2~) 
P7  hJu+Cr)2^2-  p2y2 


(12) 


where 


G(m)  = 


1-m 


(l-2m2)E*  (m)+m2H:»  (m) 


and  x  is  measured  from  the  anex  of  +h«  +,0<n 

variation  of  C  .  with  v  th©  trailing  edge.  Now  the 

of  v.  f^  P'A  7  7  «radual  for  ©  certain  range 

ation’le  negSotSf'^^ttr^?^011, ' e7^8d  *erein,  this  eeri- 
y  thSein  eat  ^1  ^  “£C'mC8l8d  18  el™>  *  »*>*“ion  (12) 


8crqmG(m)  8qmG(m) 


P,A 


(37 


(37 


(13) 


ion  (12))  by  Us  tangent 

H:e  ZTJZui^ 

SulS  SI  SU-??  adoloSea?rc“S,aSe0°n' 

ten.  ty  flow  nr  of  reference  4,  ^ch^ror^rforT4  “*  8"!°n4 


8 


flow  I:  cp  =  cx 
flow  IV :  C  = 


1 


o  <  a  <  n 
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(14) 


Thus  the  choice 


8crqmG(m) 

PV 


'IV 


8qmG(m) 


PV 


(15) 


achieves  the  desired  cancellation. 


Equations  (14)  refer  to  0°f<  J^ot^SahTof  these 
behind  the  wing  traiiinge  ge  [ ^ ce’ahead  of  the  trailing  edge  in 
2TS4JS  Til a\  <  1 f  BefSence  4  givesfor  these  disturbances 


flow  I :  C. 


,  F(ci,fc) 

I  K(k) 


flow  IV:  C  - 


|F(ri.k)-E(#,kI| 
[_  K(k)-E(h)  J 


n<  |a|  £ 


<  1 


(16) 


pitching  moment  of  the  wing. 

a.  pitching-moment  correction  ref erred  to  the  aper^ 

i.  desired,  but  it  is  more  <»£“}«*£  of AgTlh  I»ie  prooednre 
to  the  apex  of  the  t»ili«8  *ge  inasmuch  as  the  lift  is 

entails  obtaining  both  ^“^^noe  axis  back  to  the 

involved  in  the  transfer  of  the  m  jj  ls  Qf  the 

SS  TX~  C^(o)*rSThe  iXand  thi  moment  iBg»arted  by  flow  I  to 

the  elementary  hatched  trfcmgle  of  figure  3  are  therefor. 
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t 


CpPb  f(a)  da 


8c‘ 


Pb 

~2a 


iCjfCa)  i 


di 

P 


^ o 


cIP2b3f(a)  da 
24a3 

The  lift  and  the  moment  experienced  "by  the  finite  triangular  area  of 
figure  3  are  obtained  by  integrating  between  the  limits  n  and  1. 
The  function  f(a)  therein  is  the  coefficient  of  C^.  in  equa¬ 
tions  (16) «  The  substitutions  of  appendix  B  are  helpful  in  effect¬ 
ing  the  integrations.  The  results  are 

*il  r  »/s  ~i 

I  pv2  80  L  K,(n)_] 

aiM  .  _  ojsV  r  _  ,.(n)-| 

.  |  PV2  48n2  L 

The  second  of  equations  (16)  (flow  TV)  is  of  the  form 
Cp  =  Cjy-xf(a).  The  lift  and  the  moment  imparted  by  flow  IV  to  the 

elementary  hatched  triangle  of  figure  3  are  therefore 


(17) 

(18) 
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Pb 

r2  a 

CjyXf(o)  X  dx  ^ 


'-'0 


c^pVfCa)  da 


24a 


d 


Pb 

r2° 

xCjyxf (a)  x  dx 


'-'0 


da 

64  a3 


The  lift  and  the  moment  experienced  by  the  finite  triangular  area  of 
figure  3  are  obtained,  ae  before,  by  integrating  between  the  limits  n 
and  1.  The  function  f(a)  is  the  expression  in  brackets  in  the 
second  of  equations  (16).  The  result  of  the  integrations  is 


w  _ 

CjypV 

1 

E'(n)-n2K*  (n) 

|pv2 

48n2 

JL 

K'  (n)-E'(n) 

AyvM 

CjyP^b 

it  1-n2 

1  2 
o  PV 

192n3 

-L  •“  - 

4  K*  (n)-E'  (n] 

(19) 


(20) 


All  the  data  are  now  at  hand  for  calculating  the  incremental 
pitching  moment  referred  to  the  wing  apex.  First  the  incremental  lift 
will  be  needed. 


* 


c 

c 

t* 


t 


AiL  AivL 
|  pv2  |  pv2 


(21) 
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* 


Then  the  incremental  pitching  moment  is 

+  iH!!  .  Cr  _4i_  (72) 

|  PV2  |  CV2  |  PV2  |  PV2 

The  several  moments  and  lifts  in  these  two  equations  may  be 
taken  from  equations  (17)  to  (20),  using  for  Cj  and  Cjy  therein 

the  values  determined  in  equations  (15).  Then  the  left  sides  may  be 
converted  to  coefficient  form  by  dividing  equation  (21)  by  S 
and  qc/2V  and  equation  (22)  by  Sc  and  qc/2V.  Multiplication  by 
a  factor  of  2  accounts  for  both  wing  panels  on  the  complete  wing. 
The  final  results  are,  with  both  sides  multiplied  by  3, 


A(3i»Cr(m) 

r4i 

x  «/2 

2  c  A3 

4*  ^  “ 

\ 

n 

c 

K»  (n) 

3  c  n 

A3mG(m) 

»r2 

4  - 

1  */2 

n 

p2 

K'(n)_ 

E'tnj-n^^n) 
K'  (n)-E'  (n) 


(23) 


2  jr 

3  c 


£  M 

c  n 


E* (n)  E,(n)-n2K,(n) 

K'(n)  "  K'(n)-E'(n) 


2  2 
A  3 


1 


2 

«  1-n 

4  K’  (n)-E'  (n) 


(24) 


Equation  (24)  gives  the  value  of  a  correction  for  the  damping- 
in-pitch  derivative  Cm  when  the  trailing  edge  of  the  wing  is  sub- 

<1 


sonic.  This  correction  is  algebraically  additive  to  the  value 
of  CL.  computed  by  integration  of  the  basic  delta-wing  lift  distri- 
<1 

bution  alone.  Equation  (23)  gives  a  corresponding  correction  for  the 
less -important  lift-due -to -pitching  derivative  CT  .  Both  equations 

q 

are  subject  to  the  Mach  number-plan  form  limitations  expressed  by 
equations  (11a)  and  (lib). 
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Equations  (23)  and  (24)  refer  to  a  ving  pitching  about  its 
apex  with  moments  referred  to  the  apex.  The  following  transformations 
will  provide  values  of  the  corrections  to  the  derivatives  (designated 
by  primes)  referred  to  an  axis  a  distance  h  downstream  of  the  apex: 


K')  ■  *K)  ■  ®  4  W 
K) = 4K)  •  °  [Kk0-4K)] 


2h 


(25) 


The  quantities  A 


C\)  4(*\) 


on  the  right  side  are  given  by 


equations  (23)  and  (24).  The  quantities  A^BC^  ^  and  A^BCm  j  may 

be  obtained  from  reference  1  (corrected  in  the  errata  sheet  thereto). 

Similar  procedures  for  the  treatment  of  the  effect  of  subsonic 
trailing  edges  on  the  derivatives  C,  (effective  dihedral)  and  Cj 


are  considered  briefly  in  appendix  C 


P 


DISCUSSION 

The  circumstances  under  which  the  approximations  used  herein  lead 
to  acceptably  small  error  in  the  calculation  of  the  wing  loading 
require  examination.  There  are  two  points  of  departure  from  "exact" 
linearized  potential -flow  theory  employing  the  Kutta  condition  at  the 
trailing  edge:  (1)  the  basic  delta -wing  load  distribution  behind  the 
wing  is  canceled  only  approximately  by  superposition  in  each  case  of 
a  special  flow,  and  there  remains  a  generally  small  residual  lift  in 
the  wake  (fig.  2);  (2)  the  mutual  interaction  between  the  tips  and 
the  trailing  edge  (black  regions  of  fig.  1)  is  neglected. 

In  more  elaborate  treatments  (reference  1  for  angle  of  attack, 
reference  2  for  rolling  motion)  the  lift  in  the  wake  is  canceled 
almost  completely,  and  the  mutual  interaction  between  the  tip  and  the 
trailing  edge  is  also  taken  into  account.  In  reference  1  two  examples 
are  presented  from  which  the  separate  contributions  of  the  several 
degrees  of  approximation  may  be  assessed.  (Similar  examples  in  refer¬ 
ence  2  are  not  so  well  Butted  to  such  an  assessment,  because  of  an 
imperfect  correspondence.  One  of  these  examples  is,  however,  here¬ 
inafter  considered.)  The  examples  refer  to  an  untapered  wing  of 
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aspect  ratio  1.72  and  a  tapered  ving  of  aspect  ratio  3.85  and  taper 
ratio  0.179,  "both  with  63°  leading-edge  sweepback.  The  Mach  number 
is  1.5.  The  breakdown  for  lift  is  as  follows: 


Source  of  Lift 

Un tapered 
wing 

(percent) 

Tapered 

wing 

(percent) 

(l)  Basic  delta -wing  loading 

100.0 

100.0 

(2)  Tip  effect 

-14.7 

-1.6 

(3)  Symmetrical  wake  correction 

-12.1 

-2.7 

(4)  Oblique  elements  in  wake 

-.82 

-.59 

(5)  Tip-trailing-edge  mutual  interference  3.0 

•04 

Total  lift 

75.4 

95.1 

The  "symmetrical  wake  correction"  is  cancellation  flow  I  of  this 
report  (fig.  2).  The  "oblique  elements  in  the  wake"  is  an  addi¬ 
tional  flow  that,  added  to  flow  I,  (almost  exactly)  completes  the 
cancellation  of  the  basic  delta-wing  lift  behind  the  wing.  Thus  for 
the  untapered  wing  the  item  of  -0.82  percent  represents  the  error 
in  computed  wing  loading  associated  with  the  failure  of  flow  I  to 
cancel  completely  the  basic  delta -wing  loading  behind  the  wing. 
Finally,  the  item  3.0  percent  represents  the  error  associated  with 
neglect  of  the  mutual  interference  between  the  tips  and  the  trailing 
edge.  The  errors  of  the  present  treatment  therefore  amount  to 
-0.82  +  3.0  or  2.2  percent  for  the  calculation  of  the  lift  of  the 
untapered  wing.  (Note  that  these  errors  are  expressed  as  percentages 
of  item  1. ) 


The  examples  afford  a  similar  breakdown  for  the  pitching  moment 
about  the  wing  apex: 


Source  of  Moment 

Un tapered 
wing 

(percent) 

Tapered 

wing 

(percent) 

(1)  Basic  delta -wing  loading 

100.0 

(2)  Tip  effect 

-23.9 

(3)  Symmetrical  wake  correction 

-17.7 

(4)  Oblique  elements  in  wake 

-1.37 

(5)  Tip-trailing -edge  mutual  Interference 

5.3 

Total  pitching  moment 

62.3 

,  92.3 
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The  errors  of  the  present  treatment  for  the  calculation  of  the  pitch¬ 
ing  moment  for  the  untapered  wing  therefore  come  to  -1.37  +  5.3  or 
3.8  percent  and  for  the  tapered  wing  to  -0.95  percent  (again  based  on 
item  (1)  =  100  percent). 

It  is  possible  to  infer,  in  a  general  way,  from  these  examples 
the  circumstances  under  which  the  approximations  of  this  report, 
namely,  the  neglect  of  items  (4)  and  (5),  are  reasonably  satisfactory. 
The  difference  between  cancellation  flow  I  and  the  basic  delta -wing 
lift  in  figure  2  corresponds  to  item  (4).  For  a  given  section  A -A 
across  the  wake,  this  difference  decreases  with  the  ratio  Q R/PS. 

For  the  evaluation  of  item  (4)  section  A-A  may  be  taken  to  pass  through 
the  trailing  ends  of  the  wing  tips  (that  is,  QF  equals  wing  span  b 
and  P  and  S  lie  on  the  leading  edges  extended).  The  ratio  QR/PS 
for  this  section  is  designated  T.  In  the  untapered-wing 
example,  T  =  0.629;  in  the  tapered -wing  example,  T  =  0.926. 

Now  consider  variations  of  wing  sweep,  taper,  and  Mach  number 
that  leave  T  constant.  It  is  expected  that  item  (4)  will  change 
approximately  in  proportion  to  item  (3)  (exactly,  if  the  trailing - 
edge -sweep  -  Mach  number  parameter  n  is  unchanged).  On  this  basis 
item  (4)  in  the  examples  may  be  adjusted  to  equal  1  percent  and 
item  (3)  scaled  up  accordingly.  The  results  are 


Item  (3)  in 

T 

CL 

a 

a 

(percent) 

(percent) 

0*629 

-14.7 

-12.9 

•  926 

-4.6 

-4.2 

These  four  points  are  plotted  in  figure  4.  A  straight-line  variation 
of  T  against  item  (3)  is  not  inconsistent  with  the  points.  The  line 
drawn  in  is  considered  to  represent  approximately  the  condition  that 
item  (4)  equals  1  percent.  Accordingly,  for  values  of  T  below  the 
line  the  neglect  of  item  (4)  will  yield  an  error  of  less  than  1  per¬ 
cent  (based  on  basic  delta -wing  loading  =  100  percent).  The  par¬ 
ameter  T  may  be  calculated  from  the  formula 


T 


1 

m _ 

1  4  A 

—  +  ■ 

m  AB(1+X) 


(26) 
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A  corresponding  criterion  for  the  neglect  of  item  (5),  the  mutual 
interference  between  tips  and  trailing  edge *  cannot  be  formulated  with 
much  confidence.  The  ratio  of  the  area  moment  of  the  wing  regions 
affected  by  this  interference  (black  regions  in  fig.  1)  to  the  area 
moment  of  the  entire  wing  surface  is  presumed  to  be  a  rough  guide  in 
the  estimation  of  moments.  In  the  tapered -wing  example  both  item  (5) 
and  this  ratio  for  the  area  pitching  moments  are  negligible.  In  the 
un tapered -wing  example  this  ratio  for  the  area  pitching  moments  is 
17  percent  and  item  (5)  is  5.3  percent.  For  the  same  wing  this  ratio 
for  the  area  rolling  moments  is  15.4  percent  and  item  (5)  for  rolling 
motion  is  4.3  percent.  (See  subsequent  example*  taken  from  refer¬ 
ence  2.)  If  the  correspondence  can  be  stretched  somewhat*  the  rule 
that  the  percent  error  in  calculated  moment  associated  with  neglect  of 
item  (5)  is  about  one-third  the  percent  area -moment  ratio  may  be  used 
as  a  tentative  criterion.  Equations  (11a)  and  (lib)  constitute  an 
additional  limitation  that  states*  in  effect*  that  the  black  inter¬ 
ference  regions  in  figure  1  may  not  impinge  on  portions  of  the  leading 
edge.  This  limitation  was  introduced  partly  for  convenience  in  the 
integrations;  both  this  and  the  area -moment  limitation  act  in  the  same 
direction.  Of  the  two*  that  one  governs  that  results  in  the  smaller 
area -moment  ratio. 

These  several  criterions  for  the  suitability  of  the  approximations 
of  this  report  were  inferred  from  examples  worked  out  for  a  wing  at  an 
angle  of  attack.  The  similarity  of  the  features  of  the  load  cancel¬ 
lation  for  rolling  motion  and  for  pitching  motion  (fig.  2)  suggests 
that  these  criterions  are  of  the  right  order  for  these  more  general 
motions.  (Note  that  the  appropriate  axis*  in  each  case*  should  be  used 
for  the  area -moment  criterion.)  An  indication  that  the  criterion 
for  T  is  not  unconservative  for  the  case  of  rolling  comes  from  the 
following  calculation  of  Cj  for  the  un tapered  wing  of  reference  1* 

taken  (except  for  item  (3))  from  reference  2: 

Source  of  Moment 


(1)  Basic  delta -wing  loading  (reference  2) 

‘  (2)  Tip  effect  (reference  2) 

(3)  Partial  trailing-edge  correction  (equation  (9)  herein) 

(3a)  Total  trailing-edge  correction  (reference  2) 

(5)  Tip -trailing  edge  mutual  interference 

Item  (3)  is  the  approximate  trailing-edge  correction  computed  by 
equation  (9)  of  the  present  paper.  Item  (3a)  is  the  trailing-edge 
correction  computed  in  reference  2  by  substantially  exact  cancellation 


Untapered 

wing 

(percent) 

100.00 

-37.54 

-8.06 

-8.23 

4.26 
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of  the  residual  lift  in  the  wake.  Thus  the  difference  between 
items  (3a)  and  (3)  corresponds  to  item  (4)  of  the  earlier  examples. 
This  difference  is  very  small:  only  0.17  percent  of  the  uncorrected 
value  of  C,  ,  item  (1).  This  value  for  item  (4)  is  considered  to 
P 

have  some  uncertainty  in  view  of  the  fact  that  item  (3a)  was  obtained 
by  numerical  integration.  If  the  point  T  =  0.629, 
item  (3)  =  8.06  percent,  is  plotted  on  figure  4  it  falls,  as  it 
should,  within  the  region  for  which  item  (4)  is  less  than  1  percent. 

So  far  the  approximations  in  this  report  have  been  justified 
solely  on  the  basis  of  their  closeness  to  the  values  required  by  the 
"exact"  linearized  potential  theory.  For  real  wings,  however,  the 
action  of  viscosity  in  producing  a  boundary  layer  partly  vitiates 
the  predictions  of  potential  theory.  The  boundary-layer  flow  for 
sweptback  wings  is  such  that  the  deviation  from  potential  theory  is 
expected  to  be  appreciable  in  the  vicinity  of  the  wing  tips  and  the 
trailing  edge.  Thus,  the  approximations  that  have  been  introduced 
herein  are  probably  within  the  accuracy  of  applicability  of  the 
potential  theory  itself,  so  long  as  the  several  criterions  have  been 
met. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  January  6,  1950. 
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E(0,k) 

E(k) 

E*  (a)  - 


incomplete  elliptic  integral  of  second  kind  with 
amplitude  $  and  modulus  k 

complete  elliptic  integral  of  second  kind  with 


E'(n)  = 

F($,k)  incomplete  elliptic  integral  of  first  kind  with 
amplitude  0  and  modulus  k 

f(c)  function  of  a  defined  in  text 


G(m) 

h 


(l-2m2)E'(m)+m2K'  (m) 


distance  of  center  of  gravity  downstream  of  wing  apex 
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I(m) 

2(l-m^ 

(2-m2)E»(m)- 

K(k) 

complete  ell: 

modulus  k 

K'  (m)  = 

K  ( L-m2  ) 

K'  (n)  * 

k(V^F) 

P 

<L 

r 

S 

u,v,w 

V 


[■  fH 


modulus  of  elliptic  integral 
lift 


l\j  1-n2  ) 


rolling  moment  (positive  in  sense  of  right-hand  screw 
proceeding  upstream) 

pitching  moment  (positive  in  sense  of  right-hand  screw 
proceeding  spahwise  to  right) 

0  times  contangent  of  leading-edge  sweep  angle  measured 
from  y-axis 

0  times  contangent  of  trailing-edge  sweep  angle  measured 
from  y-axis 

angular  velocity  of  rolling 
angular  velocity  of  pitching 
angular  velocity  of  yawing 
wing  area 

disturbance  velocity  components  along  x-,  y- ,  «r>a  z -axe s , 
respectively 

stream  velocity 

Cartesian  coordinates:  x-axis  parallel  to  free -stream 
direction;  y-axis  horizontal  and  toward  the  right, 
looking  upstream;  z-axis  vertically  upward 


IV 


flov  IT 


1300 


1300 


NACA  TN  2146 


21 


22 


MCA  ITT  2146 


APPENDIX  C 


EXTENSION  TO  C,  AND  C, 

*3  r 

The  antisymmetric  component  of  the  basic  delta-wing  loading 
is  qualitatively  similar  for  rolling  p,  sideslip  3,  and 
yawing  r.  Thus  the  procedure  developed  herein  for  the  damping- 
in  -roll  derivative  C7  may  be  adapted  without  difficulty  to  the 

effective-dihedral  derivative  C7  and,  once  the  correct  delta- 

tp 

wing  loading  is  known  for  yawing,  to  the  derivative  C7  .  This  pro- 

*r 

cedure  will,  however,  be  less  accurate  for  C7  than  for  C7  , 

43  lV 

because  the  delta -wing  loading  to  be  canceled  in  the  wake  is  less 
well  approximated  by  flow  III.  An  approximation  that  is  closer  for 
sideslip,  but  still  simple,  is  a  flow  specified  by  u  =  Cy/(cr+x) 
in  the  region  0<la|<  n  and  w  =  0  in  the  regions  n<|c|<l. 

This  flow  is  difficult  to  solve  for  u  in  the  w  =  0  regions,  but 
a  flow  closely  similar  in  the  right-hand  w  =  0  region  may  be 
readily  solved  by  a  procedure  in  reference  4.  This  similar  flow  is 
obtained  by  replacing  w  =  0  by  u  =  0  in  the  left-hand 
region  -1  <  0<  -n .  The  final  solution  for  u  in  the  right-hand 
region  n<<7  <1  is  to  be  applied  for  both  n5L<7  <1  and  -1<  a  <-r_. 
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Figure  2.  -  Sweptback  wing  showing  basic  delta-wing  load 

distribution  along  section  A-A  and  approximate  cancellation 
of  this  load  behind  wing  by  superposition  of  special  flows. 
Special  flows  are  plotted  with  reversed  sign.  Cross- 
hatched  areas  represent  induced  changes  in  loading  ahead 
of  trailing  edge. 
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Trai ling-edge  correction,  item  3,  percent 


Figure  4.  -  Estimated  variation  of  t  with  trailing-edge  correction 
computed  by  present  method  under  condition  that  error  in  correc¬ 
tion,  exclusive  of  tip  interaction,  equals  1  percent  of  contribu¬ 
tion  of  basic  delta-wing  loading. 
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